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1 Introduction 

The theory of wave scattering by small bodies was initiated by Rayleigh (1871). Thomp- 
son (1893) was the first to understand the role of magnetic dipole radiation. Since then, 
many papers have been published on the subject because of its importance in applica- 
tions. Prom a theoretical point of view there are two directions of investigation: 

(i) to prove that the scattering amplitude can be expanded in powers of ka, where 

k = 2it/X and a is a characteristic dimension of a small body, 

(ii) to find the coefficients of the expansion efficiently. 

Stevenson (1953), Senior and Kleinman can be mentioned among contributors to the 
first topic. 

To my knowledge, there were no results concerning the second topic for bodies of an 
arbitrary shape. Such results are of interest in geophysics, radiophysics, optics, colloidal 
chemistry and solid state theory. 

In this paper we review the results of the author on the theory of scalar and vector 
wave scattering by small bodies of an arbitrary shape with the emphasis on practical 
applicability of the formulas obtained and on the mathematical rigor of the theory. For 
the scalar wave scattering by a single body, our main results can be described as follows: 

(1) Analytical formulas for the scattering amplitude for a small body of an arbitrary 
shape are obtained; dependence of the scattering amplitude on the boundary con- 
ditions is described. 
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(2) An analytical formula for the scattering matrix for electromagnetic wave scattering 

by a small body of an arbitrary shape is given. Applications of these results are out- 
lined (calculation of the properties of a rarefied medium; inverse radio measurement 
problem; formulas for the polarization tensors and capacitance). 

(3) The multi-particle scattering problem is analyzed and interaction of the scattered 

waves is taken into account. For the self-consistent field in a medium consisting of 
many particles (~ 10^^), integral-difi'erential equations are found. The equations 
depend on the boundary conditions on the particle surfaces. These equations offer 
a possibility of solving the inverse problem of finding the medium properties from 
the scattering data. For about 5 to 10 bodies the fundamental integral equations of 
the theory can be solved numerically to study the interaction between the bodies. 

In section 2 the results concerning the scalar wave scattering are described. In section 
3 the electromagnetic scattering is studied and the solution of the inverse problem of radio 
measurements is outlined. In section 4 the many-body problem is examined. 

2 Scalar Wave Scattering by a Single Body 

Consider the problem 

'duo 



dN 



hu( 



exp (i/c a; ) ^/ , x , , x , . 

V ~ j— — -jin^k) as \x\ — > oo, — = n, (2) 
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where Q = R^\D, D is a bounded domain with a smooth boundary F, N is the outer 
normal to F, Uq is the initial field which is usually taken in the form Uq — exp {ik{i', x)}. 
We look for a solution of the problem (l)-(2) of the form 

exp (/7,7- ./) a(f ) , , 

and for the scattering amplitude / we have the formula 

f = — ex.p{-ik{n, t)}a{t, k)dt = — ao{t)dt + 0{ka), (4) 

4:71 Jy 4:71 Jy 



where 



a{t, k) = aoit) + ika^it) + ^^(J2{t) + ... (5) 
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Putting (3) in the boundary conditions (1) we get the integral equation for cr: 

.duo 



a = A{k)a - hT{k)a - 2huo + 2 



ON' 



where 



A{k)a 



(6) 



(7) 



Expanding a, A{k) and T(/c) in the powers k and equating the corresponding terms in (6) 
we obtain, for h = 0, i.e. for the Neumann boundary condition, the following equations: 

CTq = AqCTo, (8) 



(Ti = AqCTi + Aiao + 2 
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0-2 = ^0(72 + 2Ai(Ti + ^2(70 + 2' 
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etc, where 

(ik)'^ {iky 
A{k) = Ao + ikAi H + • • • ; mq = moo + ^^mqi H ^"02 + • • 

Expanding / in formula (4) we obtain, up to the terms of the second order: 
/ = h !t ^odt + ik{^ /p aidt + (n, /p aoit)tdt) } 



(ik)^ r 1 



{ 4^ /p cr2dt + ^ (n, /p aitrft) + ^ /p (To(n, 



(10) 



'12) 



From (8) it follows that ctq = and from (9) it follows that Jp aidt = 0. Some calculations 
lead to the following final result (see [|17|): 



ikV 
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f 4^ ^Pi^'Pdx^ 
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(13) 



z=0 



Usually Mo = exp{ifc(z/, z)}, and in this case formula (13) can be written as: 

k^V 



f 



{PpqVqHp + 1) 



:i3' 



where over the repeated indices the summation is understood, V is the volume of the 
body D and j3pq is the magnetic polarizability tensor of D. Note that / ~ k'^a? the 
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scattering is anisotropic and is defined by the tensor j3pq. Formula (25) below allows one 
to calculate (3pq. 

For h = oo (the Dirichet boundary condition) integral equation (6) takes the form: 

T{k)a = -2uq. (14) 

Hence 

a^dt 



An 



-no 



St 



Since ka << 1 the field uo |r = uq{x, /c)|^=o' where the origin is assumed to be inside the 
body D. From the above equation it follows, that J-^aodt = —Cuo, 

where C is the capacitance of a conductor with the shape D. Hence for the Dirichet 
boundary condition, f ^ a, where a is a characteristic length of D, and the scattering is 
isotropic. 

For h 0, using the same hne of arguments, it is possible to obtain the following 
approximate formula for the scattering amplitude: 

^"-MTTXscm) (^6) 

where S = meas (F), i.e., the surface area of F, and C is the capacitance oi D. U h is 
very small {h ~ k'^a^) the formula for / should be changed and the terms analogous to 
(13) should be taken into account. 



3 Electromagnetic Wave Scattering by a Single Body 



If a homogeneous body D with the parameters e, /i, a, is placed into a homogeneous 
medium with the parameters Eq, /iq, c"o; then the following formula for the scattering 
matrix S was established by the author (see ||T7[| ): 



S 



k^V 
Air 



/^o Ai + 0^22 COS 9 — sin 9, 021 cos 9 — a^i sin 9 — /io/5i2 
"12 - /^o/32i COS 9 + Ho/3u sin 9, an + /io/?22 cos 9 - /io/532 sin 9 



where S is defined by the formula 

= S 



f2 
fl 
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(17) 



9 is the angle of scattering, Ei, E2 are the components of the initial field, /i, /2 are the 
components of the scattered field in the far field region multiplied by exp(i/c|a;|). 
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the plane YOZ is the plane of scattering, aij — aij('y), is the polarizabihty tensor, 
^ = {e — Sq) / {e + So) and f3ij = 1) is the magnetic polarizabihty tensor. 

If one knows S one can find all the values of interest to physicists for electromagnetic 
wave propagation in a rarefied medium consisting of small bodies. The tensor of refraction 
coefficient can be calculated by the formula riij — 6ij + 27rNk~'^Sij{0), where N is the 
number of bodies per unit volume. The tensor Oiij{'y) can be calculated analytically by 
the formula 



,(7) 



(")/ 



(7) <Aq'\ 0<q<l, 



(19) 



where A, and q are some constants depending only on the geometry of the surface, and 
«S?-^Efe^ ^ ~] ht\ n>l. (20) 

In (20) 



m=0 



(1) 




N,{s)Nj{t)dsdt 



(21) 



r Jr 



= Ir Ir dsdtNi{s)Nj{t) J ' ' ' J TTt"^ (^i' ^) • • • ^ (^"^-i' ^"^-2) dty . . . 

(22) 



m— 1 



dti . . . dt^_i; ^(t,s) = g|- ^. 



In particular 



2,(1) 1^(1) 

4'*W=2(7 + T^)*«-^, = (23) 



For particles with 11 = 110 and e not very large, so that the depth 6 of the skin layer is 
considerably larger than a, one can neglect the magnetic dipole radiation and in formula 
(17) for the scattering matrix one can omit the terms with the multipliers 

The vectors of electric P and magnetic M polarizations can be found by the following 
formulas, respectively, 

P, = a,,{^)VeoE,, 7:=^, (24) 
where Ej is the initial field, over the repeated indices one sums up, and 

Mi = aij{^)VnQHj, +pijViJ,oHj, 7 := ^ := aij{-l), (25) 
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where Hj is the initial field and the second term on the right hand side of equality (25) 
should be omitted if the skin- layer depth 5 » a. 

The scattering amplitudes can be found from the formulae 

fE = ^[nAPM] + ^\[f[MM. (26) 
47r£o 47r V Sq 



fH=J-[n,fE], (27) 




where [A, B] stands for the vector product A x B, and P and M can be calculated by 
formulae (24), (25), (19), (20), (22). If 5 » a one can neglect the second term on the 
right-hand side of (26). 

It is possible to give a simple solution to the following inverse problem which can be 
called the inverse problem of radiomeasurements. 

Suppose an initial electromagnetic field is scattered by a small probe. Assume that the 
scattered field E', H' can be measured in the far field region. The problem is: calculate 
the initial field at the point where the small probe detects H' . This problem is of 
interest, for example, when one wants to determine the electromagnetic field distribution 
in an antenna's aperture. Let us assume for simplicity that for the probe 5 >> a, so that 

r 47r£o 

From (28) one can find P — Ui (P, rii) = E' (ni) b where b — ^^^^^^''^ A measure- 
ment in the n2 direction, where (711,^2) = 0, results in P — n2(P, ^2) = £"(^2)6. Hence 
(ni,P) = b{E'{n2),ni). Thus P = b{E' (rii) + ni{E' (ui) , rii)} . But 

Pi = aij{-f)VeoEj. (*) 

Since V and Sq are known and aij{j) can be calculated by formulae (19), (20) and the 
matrix cty is positive definite (because ^aijVeoEjEi is the energy) it follows that system 
(*) is uniquely solvable. Its solution is the desired vector E. 

Let us give a formula for the capacitance of a conductor D of an arbitrary shape, 
which proved to be very useful in practice: 



C(") = AttsS' I 1^ JJ^^J^...J^ ^(t, h)... i;{tn-i, tn)dh ...dtA (29) 

V, n times / 

5 

C(°) = ^<C, J^ff^, 5 = mea.F. (30) 
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It can be proved that 

< Ag", 0<g<l, (30') 
where A and q are constants which depend only on the geometry of F. 
Remark 3.1. The theory is also applicable for small layered bodies. 

Remark 3.2. Two sided variational estimates for aij and C were given in and /jl^. 

4 Many Body Wave Scattering 

First we describe a method for solving the scattering problem for r bodies, r ~ 5 — 10, and 
then we derive an integral-differential equation for the self-consistent field in a medium 
consisting of many (r ~ 10^^) small bodies. We look for a solution of the scalar wave 
scattering problem of the form 

u = uo + y f ^-^^^^a,{t)dt. (31) 



Applying the boundary condition, 



u 



1 < J < r, (32) 



we obtain the system of r integral equations for the r unknown functions aj. In general 
this system can be solved numerically. When d << X, where d = mini-^jdij, and dij is 
the distance between i — th and j — th body, the system of the integral equations has 
dominant diagonal terms and it can be easily solved by an iterative process, the zero 
approximation being the initial field uq- 

U ka >> 1, d » a, but not necessarily d » X, the average field in the medium 
consisting of small particles can be found from the integral equation 

u{x,k) = Uo{x)- [ ^^^^l^—^q{y)u{y,k)dy. (33) 

Here q{y) is the average value of hjSj (l + hjSjCj^^ ^ over the volume dy in a neigh- 
borhood of y for bodies with impedance boundary conditions. For hj = oo (the Dirichlet 
boundary condition) and identical bodies, one has q{y) = N{y)C, where N{y) is the num- 
ber of the bodies per unit volume and C is the capacitance of a body. For the Neumann 
boundary condition the corresponding equation is the integral-differential equation: 

u{x,k)=uo{x,k)+ [ ^J^^^fs^y^du^ ^^^^ + Uu{y,k)b{y)]dy, 

(34) 
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where 



h{y)=N{y)V, B„iy) = tkVP,,N{y), (35) 

V is the volume of a body, and f3pg is its magnetic polarizabihty tensor (see forumla 
(25)). The solution to equations (33) and (34) can be considered as the self-consistent 
(effective) field acting in the medium. 

Equations (33) and (34) allow one to solve the inverse problems of the determination 
of the medium properties from the scattering data. For example, from (33) it follows 
that the scattering amplitude has the form 

For a rarefied medium it is reasonable to replace u by Uq (the Born approximation) and 
to obtain 

f^^-rl exp{-iA;(n,y)}?(y)xio(?/,A;)(iy. (37) 

47r Jr3 

If uq — exp{i/i;(i/, x)} formula (37) is valid for k » 1 with the error 0{k~^) if 

|g(a;)| + |Vg(x)|<^^-^p^, £>0. 

Hence if / is known for all < A; < cxo, and all {n — v) E S2, where 5*2 is the unit 
sphere in R^, the Fourier transform of q{y) is known, and q can be uniquely determined. 
If q{y) is compactly supported, i.e. q{y) = outside some bounded domain, then / 
is an entire function of k, and knowing / in any interval [ko,ki], < ko < ki, for all 
{n — 1/) e S2 one can find / for all < /c < 00 uniquely by analytic continuation, and 
thus one can determine q{y) uniquely. 

Let us consider the r-body problem for a few bodies, r ~ 10 (small r). Assume that 
the Dirichlet boundary condition holds. Let us look for a solution of the form 

u{x) =uo + y f ^^P(^^I^-^I) (^)^^. (38) 
Jr 47r b — t\ 

The scattering amplitude is equal to 

1 f 

f{n,k) = — y^exp{—ik(n,tj)} / ex.p {—k(n,t — tj)} aj{t)dt, (39) 

where tj is some point inside the j-th body. Since ka « 1 this formula can be rewritten 
as: 

1 

f{n, k) = —J2 {-ik{n, tj)} Qj, (40) 
^ 3=1 
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where 

' • (41) 



Qj = I ajodt + 0{ka), ajo = (Jj 



k=0 



This is the same hne of arguments as in Section 2. Using the boundary condition one 
gets: 

f exp{ik\xm - 1\) f exp{ik\xn, - 1\) 

With the accuracy of 0{ka), this can be written as: 

where cimj = l^^m — If C*™ is the capacitance of the m-th body one can rewrite (43) 
as: 

This is a hnear system from which Qj, 1 < m < r, can be found. If dmjC^ » 1, 
this system can be easily solved by an iterative process. If {Qj} are known, then the 
scattering amplitude can be found from (40). 

More details about the described theory the reader can find in the References, espe- 
cially in he monograph |[T7[|. 
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